Topological Flat Bands from Dipolar Spin Systems 
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We propose and analyze a physical system that naturally admits two-dimensional topological 
nearly flat bands. Our approach utilizes an array of three- level dipoles (effective S = 1 spins) 
driven by inhomogeneous electromagnetic fields. The dipolar interactions produce arbitrary uniform 
background gauge fields for an effective collection of conserved hardcore bosons, namely, the dressed 
spin-flips. These gauge fields result in topological band structures, whose bandgap can be larger 
than the corresponding bandwidth. Exact diagonalization of the full interacting Hamiltonian at 
half-filling reveals the existence of superfluid, crystalline, and supersolid phases. An experimental 
realization using either ultra-cold polar molecules or spins in the solid state is considered. 



Single-particle flat bands, where kinetic energy is 
quenched relative to the scale of interactions, are being 
actively explored in the quest for novel strongly corre- 
lated phases of matter [1, 2]. Prompted by the analogy 
to Landau levels, recent efforts have focused on topo- 
logical flat bands (TFB) - lattice models in which the 
band-structure also harbors a non-trivial Chern invari- 
ant; interacting particles living in such a kinetic environ- 
ment are leading candidates for the realization of frac- 
tional Chern insulators, the lattice version of fractional 
Hall states [1-8]. In addition to their fundamental sig- 
nificance, such fractional topological insulators may also 
enable the realization of robust quantum memories [9], 
protected state transfer [10], and topological quantum 
computing [11]. Seminal recent work has highlighted that 
certain classes of highly-engineered two-dimensional tight 
binding models can indeed exhibit topological nearly flat 
bands [12-16]. However, the identification of a physical 
system whose microscopies naturally admit TFB remains 
an outstanding challenge. 

In this Letter, we demonstrate the emergence of syn- 
thetic gauge fields for an ensemble of interacting hard- 
core bosons — the effective spin-flips of pinned, three- 
level dipoles in a two-dimensional lattice. Underlying 
these gauge fields are two key ingredients: spatially vary- 
ing, elliptically-polarized external (microwave or optical) 
fields break time-reversal symmetry, while anisotropic 
dipolar interactions induce orientation-dependent phases 
onto the hopping hardcore bosons. The combination of 
these effects naturally produces nontrivial Chern num- 
bers in the band structure and, when tuned appropri- 
ately, results in the emergence of flat bands due to hop- 
ping interference. 

Let us consider a square lattice composed of fixed, 
three-state magnetic or electric dipoles placed in a static 
external field. Such an arrangement naturally arises 
in experimental systems ranging from ultra-cold polar 
molecules [18-24] and Rydberg atoms [25-27] to solid- 




FIG. 1. Schematic representation of a 2D dipolar droplet. 
The grey droplet represents a 2D array of interacting tilted 
dipoles. The dipoles are tilted by a static field in the z direc- 
tion, oriented at Go, $o relative to the lattice basis {X, Y, Z}. 
Hij is a vector connecting dipoles in the XY plane. 



state spins [28, 29] and magnetic atoms [30]. As shown 
in Fig. 1, the dipoles occupy the {X, Y} plane and couple 
via dipole-dipole interactions, 
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where k is l/47reo for electric dipoles or /io/47r for mag- 
netic dipoles, and R^- connects the dipoles and dj. 
The three states of each dipole, which we label as |0), 
| ± 1), are eigenvectors of the ^-component of (rotational 
or spin) angular momentum. We assume that the | =b 1) 
states are degenerate while the |0) state is energetically 
separated from them (Fig. 2a). 

Each three-level dipole is driven by electromagnetic 
fields of Rabi frequency Q- and detuning A as shown 
schematically in Fig. 2a. With |fi+|, <C A, the ap- 
proximate eigenstates (dressed states) are: |0), \B) — 
a(\ - 1) + /3|1», and \D) = a*(-/3*\ - 1) + |1)), where 
a = ap = and ft = Vl^-I 2 + l^+l 2 - 

The energies of these dressed states are Eq = —Vt 2 / A, 
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FIG. 2. a) Depicts the on-site level structure and the two- 
photon driving scheme. £7_(r) is left-circularly polarized, 
while is right-circularly polarized. These levels could, 

for example, be adiabatically connected to the J — 1 mani- 
fold of a rigid rotor as one turns on a DC electric field (see 
Eq. (6)). b) Square lattice with a single tilted dipole per ver- 
tex. We index columns of the lattice by £ and plaquettes by 
pi. For a particle traversing the edge of a single plaquette, 
there are two contributions te and t' e to W(pe); each contribu- 
tion occurs twice as represented by the red and blue colored 
arrows. A simple periodic gradient of f3 enables uniform n/N 
flux per plaquette. 



E B = A + f2 2 /A, and E D = A respectively. We let 
d represent the typical size of the dipole moment and 
Ro be the nearest-neighbor spacing; by ensuring that 
nd 2 /Rq <C ft 2 / A and so long as we initially avoid pop- 
ulating \D), the system remains within the subspace lo- 
cally spanned by |0) and \B). 

Thus, it is natural to view \B) as representing an ef- 
fective hardcore bosonic excitation (spin-flip), while |0) 
represents the absence of such an excitation. Recasting 
this system in terms of operators a\ = \B)(0\i (rii = ajoj) 
yields a 2D model of conserved hardcore lattice bosons, 
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where we define the hopping Uj = — (BiQj\Hdd\OiBj), 
the on-site potential tu = ^2j^i((0i0j\H dd \0i0j) - 
(BiOj\Hdd\BiOj)) : and the interaction Vij = 
(BiBj\H dd \BiBj) + (OiOjlHddlOiOj) - {B^H^B^) - 
(0iBj\H dd \0iBj). The conservation of total boson 
number, Ni = X^ a ! a i> arises from the condition 
nd 2 /Rq <C A, which ensures that particle- number non- 
conserving terms of H dd are energetically disallowed. 

The functional form of the effective hardcore bosonic 
Hamiltonian Eq. (2) arises for any system of pinned, 
three- level dipoles. The parameters in Hb are given by 
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where d° (d B ) is the permanent i-dipole moment of the 
|0) (| B}) state, doi is the transition dipole moment from 
|1) to |0) [31], Xi ~ a i(l,fti) T is the normalized drive- 
spinor on site z, qo = ^(1 — 3cos 2 (<£ — $0) sin 2 (6o)), 
q 2 = _|[ CO s($ - $o)cos9 - isin($ - $o)] 2 , & are the 
Pauli matrices, and (R, <£) is the separation in po- 
lar coordinates (Fig. 1). We have suppressed the explicit 
ij dependence of R, <3>, qo, and q 2 . While the form of 
df, and hence of interactions, depends on the underly- 
ing implementation, the single-particle band structures 
that can be achieved via driving are independent of such 
details [32]. 

Let us first explore these topological single-particle 
bands and illustrate the interplay between the driven 
breaking of time-reversal and the anisotropic dipolar in- 
teraction. As a simple example, we demonstrate how to 
achieve a synthetic background gauge field with uniform 
flux 7t/N per plaquette on a square lattice (assuming only 
nearest-neighbor hops). We choose the "magic" electric 
field tilt, (e ,*o) = (sin^V^),?^), where q = 
along X and Y. This choice allows us to isolate the 
terms of H dd that harbor intrinsic phases, namely, those 
associated with dj d^ and d~dj [21, 31]. Moreover, it 
simplifies the form of nearest-neighbor hopping to 
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The microscopic breaking of time-reversal arises from the 
asymmetry between left- and right- circularly polarized 
radiation and is captured by the ratio f3 = In- 
deed, physically, it is this ratio that defines each hardcore 
boson | by setting the relative admixture between the 
|1) and | — 1) states. Keeping j3 real, let us now con- 
sider varying the intensities of the drive fields along the 

= 7r/4 direction in a periodic fashion. 

For each plaquette, we define the Wilson loop, W(p) = 
Yidp •> which is identical along columns indexed by £ 
(Fig. 2b). The flux in a plaquette is then the phase of this 
Wilson loop, — arg[W(p^)] = arg[tftg 2 ], where tg are t' £ 
are the hops depicted in Fig. 2b. Taking 0£ = arg(^) and 
noting that 6' £ = arg(t^) = — yields the phase of the 
Wilson loop as = 20 £ — 20£+\. To achieve a uniform 
tt/N flux per plaquette, we can take 0i+\ = r] — £^, 
where 77 £ R is a constant to be specified. From the 
definition of 0£, one finds a simple recursion relation for 
(3, 



Pt+i = sin(f + 
Pi sin(f+7?-£^) 



(5) 
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with maximum periodicity AN [48]. Starting from any 
initial j3\, Eq. (5) yields a recursively generated drive pat- 
tern which achieves the desired uniform tt/N background 
gauge field. 
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FIG. 3. a) Schematic representation of the two-site unit cell 
lattice with f3 = /3±, @2- The dotted box outlines a single unit 
cell. There is a flux <1>, — which alternates in neighboring 
square plaquettes. The direct lattice vectors g\ and are 
depicted as purple arrows. While all hops are present with 
amplitude decaying as 1/i? 3 , only nearest-neighbor (solid) 
and next-nearest-neighbor (dashed) hops are shown, b) The 
topology of bulk bands as a function of complex /?2 for /3i G M. 



While the uniform flux per plaquette is reminiscent of 
the square lattice Hofstadter problem [33], we emphasize 
that the physics of these driven dipoles is significantly 
richer, owing to the additional modulation of Uj. The 
background flux field arises, in part, from the natural 
phases associated with the dipolar interaction. This en- 
sures that (as in [34]) the number of flux quanta per 
plaquette is not limited by the magnitude of laser inten- 
sities, contrasting with the majority of previous synthetic 
gauge field proposals, where the scaling to high artificial 
fluxes is extremely difficult [35-39]. 

To illustrate the symmetry breaking required for the 
generation of gapped Chern bands, we now turn to a de- 
tailed study of Hb restricted to a two-site unit cell (re- 
maining at the "magic" tilt), as depicted in Fig. 3a. This 
restriction has the virtue of being analytically tractable 
and allows us to identify the anti-unitary symmetries as- 
sociated with the Dirac points [40, 41]. Let us consider 
^ = on the two sites of the unit cell and include 

all terms up to next-next-nearest neighbor. The topol- 
ogy of the bands depends on the relative ratio of f3\ and 
/?2- For /3i G M, the phase diagram in Fig. 3b illustrates 
the Chern invariant of the bottom band as a function of 
the complex /^-plane. There exist two circles of gapless 
(Dirac) points protected by distinct anti-unitary symme- 
tries. 

On the outer circle, the Hamiltonian harbors an effec- 
tive time-reversal symmetry evidenced by the fact that 
all Wilson loops are real. On the inner circle, the Hamil- 
tonian is invariant under conjugation combined with in- 
version through the mid-point of an edge. Away from 
the magic tilt and including longer range terms, both 
circles deform smoothly. For generic /3\ and /?2, these 
Dirac points gap in a manner that causes the bands to 
acquire topology (Fig. 3b). Underlying these non-trivial 
Chern bands is the aforementioned synthetic gauge field 
associated with the effective dressed spin- flips. 

Implementation — An experimental realization of our 



proposal can be envisioned with systems including either 
electric (e.g. polar molecules) or magnetic (e.g. solid-state 
spins) dipoles. As previously mentioned, the form of df 
depends on this choice, since the permanent dipole mo- 
ment of the | =b 1) states have either the same or opposite 
signs. We emphasize that the long intrinsic lifetimes of 
such systems make them ideal for the consideration of 
driven, non-equilibrium phenomena [42, 43]. 

To be specific, we now focus on diatomic polar 
molecules (trapped in a deep optical lattice) in their 
electronic and vibrational ground state. We utilize mi- 
crowave fields to dress the molecules and partially polar- 
ize them with an applied DC electric field along z (Fig. 1); 
ignoring electronic and nuclear spins, this yields a single- 
molecule Hamiltonian, 

H m = BJ 2 — d z E + H.MW-, (6) 

where B is the rotational constant, J is the rotational 
angular momentum operator, d z is the z component of 
the dipole operator, E is the magnitude of the applied 
DC field, and Hmw characterizes the microwave dressing 
[21,42]. 

In the absence of applied fields, each molecule pos- 
sesses rigid rotor eigenstates | J, M). The applied electric 
field E mixes eigenstates with the same M, splitting the 
degeneracy within each J manifold and inducing a finite 
permanent dipole moment for each perturbed rotational 
state. We choose from among these states to form the 
effective three- level dipole; an example of one possibility 
for |0), | d= 1) is shown in Fig. 2a. Since these | ± 1) states 
have an identical induced dipole moment d 1 , one finds 
that df = d 1 , and hence, 

ViJ =2^(d°-dr. (7) 

The relative strength of the interaction Vij/Uj is thus set 
by (d° — d 1 ) 2 /d 2 )1 ; this is a highly tunable parameter and 
can easily reach ~ 100 for certain choices of rotational 
states and DC electric field strengths [21]. This tun- 
ability suggests the possibility of observing interaction- 
dominated many-body phases in such a polar-molecule 
implementation. 

The main challenge in an experimental realization of 
our proposal lies in the spatial modulation of the drive 
fields at lattice scale. For spins in the solid-state and 
on-chip polar molecule experiments, one might envision 
using near-field techniques. A more straightforward ap- 
proach, suitable for molecules, is to utilize pairs of opti- 
cal Raman beams [44]. For example, the so-called lin±lin 
configuration [45] automatically ensures that Cl and A are 
identical on all sites and moreover, generically produces 
gapped topological band-structures. 

Many-body phases — To illustrate the power of the 
present approach, we briefly explore two examples of cor- 
related ground state phases which arise in the Hamilto- 
nian Eq. (2). As Hb conserves boson number TV, we may 
consider its many-body physics at finite filling fractions 
v (particle number per unit cell). Let us work with a 
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FIG. 4. Phase Transitions in topological flat bands of 2D 
driven dipoles. a) Band structure for (Bo,^o) = (0.46,0.42), 
p 1 = 3.6e 2 69i , and f3 2 = 5.8e 5 - 63 \ Significantly flatter 
band structures with flatness ratio > 10 can be obtained for 
slightly generalized level configurations [44] . b) Structure fac- 
tor S(R,0) = (n(R)n(0)) for filling v = 1/2 in KMS [48] and 
c) SSS regime; size of circles indicates weight, d) Spectral gap 
density plot as a function of varying MW drive for parameters: 
(B ,$o) = (0.66,7r/4), ft = -2.82e i01 , 2 = -4.84e~^ 2 and 
(d° - cZ 1 ) 2 /^! « 2.8. The transition from the SF, which 
has a unique finite-size ground state, to the degenerate SSS 
shows as a collapse of this gap. e) Spectral flow in the ground 
state momentum sector of the SSS under twisting of the bo- 
son boundary condition in the g\ and f) g 2 directions. For 
the N s — 24 lattice with 6 bosons, momentum sectors return 
to themselves after 2n in 0i and after 4tv in 2 . 



two-site unit cell and truncate the dipolar interactions 
at next-next-nearest-neighbor order. Bosons residing in 
a strongly dispersing band structure generically form su- 
perfluids in order to minimize their kinetic energy. In- 
teraction dominated phases arise when the single-particle 
bands disperse less than the scale of interactions. Numer- 
ical optimization of the flatness ratio (bandgap/lowest 
bandwidth) over the six-dimensional parameter space of 
microwave driving and tilt angle reveals approximately 
flat Chern bands in several regions of phase space. The 
flatness of these bands (Fig. 4a) derives from interfer- 
ence between the hopping in different directions and, mi- 
croscopically, owes to an interplay between the natural 
anisotropy associated with dipolar interactions and the 
spatial variation of the drive fields. These quenched band 
structures provide a fertile starting point for discovering 
strongly correlated phases [12-16]. 

As a first example, we consider the band structure 
depicted in Fig. 4a, where the lower/upper band carry 
Chern index, c = =pl (parameters in caption). Ex- 
act diagonalization at filling fraction v = 1/2 and rel- 
ative interaction strength (d° — d 1 ) 2 /d^ 1 « 6 reveals a 
knight's move solid (KMS) phase with a 4-fold degen- 



erate, gapped, ground state. The real-space structure 
factor S(R,0) = (n(R)n(0)} (at total number of sites, 
N s = 32) in Fig. 4b illustrates the knight's move rela- 
tionship of the bosons in the ground state. Twisting the 
boundary condition of the KMS in the g^ directions 
(Fig. 4b) does not significantly affect the ground state 
energy, as expected of an insulator [48]. 

Many other commensurate phases arise as as we tune 
the driving fields to other regions of phase space. Fig- 
ure 4d shows a phase diagram containing both superfluid 
(SF) and striped supersolid (SSS) phases. We can charac- 
terize the SSS arising at </>i = <p2 = 0.1 as follows: First, 
diagonalization reveals the existence of three degenerate 
ground states in the sectors: hi = 0, k\ = 0, 27r/3, 47r/3. 
Consistent with striped ordering, the structure factor 
shows density stripes in the c/2 direction (Fig. 4c). How- 
ever, each of these stripes has incommensurate boson 
number, suggesting derealization along the stripes. To 
wit, for N s — 24, the 6 hardcore bosons are distributed 
evenly along two stripes, each containing 4 sites. Strong 
phase coherence along the stripes shows up in the sensi- 
tivity to twists in the cj2 direction, while transverse twists 
produce essentially no dispersion, as shown in Fig. 4e,f. 
Enlarging the system to N s — 36 does not alter this qual- 
itative picture. We note that the stripe ordering breaks 
the underlying translational symmetry of the lattice but 
not any existing rotational symmetry, as these are explic- 
itly broken by the tilt of the dipoles. 

Conclusion — The above considerations indicate that 
driven dipolar spin arrays can effectively generate topo- 
logical nearly flat bands. Implementations using either 
diatomic polar molecules in optical lattices [18-24] or lo- 
calized spins in the solid-state (such as those associated 
with NV defect centers [28]) naturally yield examples of 
interaction-dominated phases. 

Our method also opens the door for a number of in- 
triguing directions. In particular, the adiabatic prepara- 
tion and detection of single-excitation states may provide 
an elegant approach to probing chiral dynamics, edge 
modes, and the Chern index [10, 44, 46]. More gener- 
ally, dynamical preparation, manipulation and detection 
of many-body states in such driven topological systems 
remains an exciting open question [47] . Finally, the large 
available parameter space holds the promise of more ex- 
otic phases, such as fractional Chern insulators [44]; real- 
izing such phases in an effective dipolar spin system may 
provide a building block for topological quantum infor- 
mation processing. 
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Proof of uniform n/N: Here, we prove that a periodic drive pattern enables arbitrary 
uniform tt/N flux per square plaquette, as shown in Fig. 2b of the main text. The flux in 
each plaquette p is the phase of the Wilson loop, W(p) = Yidp^j- For a given plaquette, 
pi, <^ = arg[ty(^)] = &rg[tjtf ], where t £ = fa-iw + faw, t\ = few + fy+iw (dropping the 
normalization factor which will not contribute to the argument) and w = e l7T / 3 . Taking 0£ = 
arg(i^) and & t = arg(i^) yields the phase of the Wilson loop as $^ = 20^+20^ = 20^ — 20^ +i . To 
achieve a uniform tt/N flux per plaquette, we can take 0^ + i = rj — £^ + a£7T (note that in the 
main text the a £ term was dropped for notational simplicity), where a £ G {0, 1} and rj G R 
is a constant to be specified. Defining T£ = fy+i/fy yields = Ylk=^\. r kfti, implying that the 
periodicity of (3 is L such that Yl k=1 r k = 1. We now prove, for all TV G Z >0 , that there exists 
a periodic choice of f3 with a maximum periodicity of AN which generates the desired uniform 

B 

background flux. First, note that ti+i/fy = r^w + w and hence, = arg(r£W + w) + 6^7r, 
where hi G {0, 1}. Simple geometry reveals that taking 
& _ = A±i = ^(1-77 + ^) 

i— h 'A sin(f +17-^) 1 J 

> 

ON ensures that arg(r^u> + w) = 77 — £^ + Q 7 ^ where q G {0, 1}. Thus, we have succeeded 

r-- 

in achieving 0£ of the desired form with ci£ = be + Q(mod 2). Inspection reveals that for 

^ N e 3Z, ElS gfggi = 1 for all V e K, while for iV e Z\3Z, IlS SfS^fj = 1 

^ rj = Xtt/N with integer A. □ 

As a simple example, to generate a uniform n/3 background gauge field, one can utilize a 
drive pattern with period 2N = 6; with f3 = {1, 5, —10, |, |, |} (e.g. 77 = tan _1 (^g) — 7r/6), 
one finds that the topological bandstructure exhibits six bands, where the lowest carries 
Chern number c = — 1. 



Knight's Move Solid: The knight's move solid discussed in the maintext exhibits 4 de- 
generate ground states in momentum sectors k\ = 0, k 2 = 0, 7r/2, 7t, 3tt/2. As shown in 
Fig. Sl(b), twisting the boundary condition in the g 2 directions (Fig. 4b) does not af- 
fect the ground state energy consistent with an insulator. By enlarging the system size, 
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FIG. SI: Knight's move solid, a) Exact diagonalization at N s = 24 reveals a spectral gap 
density plot, which varies as a function of MW drive amplitude for (60, $0) = (0.21, 0.42), 
pi = 2.8e L69i , and (3 2 = 2.8e 5 - 63 \ Although the MW parameters are slightly different than in the 
main text, the qualitative KMS phase is unchanged; this choice of parameters better illustrates 
the transition from a SF, which has a unique finite-size ground state, to the degenerate KMS. b) 
Spectral flow of the KMS ground state in the k\ — k 2 = momentum sector under twisting the 
boson boundary condition in the g± and g 2 directions for parameters indicated by the white circle 
in a). For the N s — 24 (3 x 4 x 2) lattice with 6 bosons, momentum sectors return to themselves 
after 2tv in 0\ and after 4tv in 9 2 • In the KMS phase, the ground state energy is insensitive to 
twisting, (c) This contrasts with the quadratic response of the energy in the phase-coherent SF. 
Depicted is the spectral flow of the k\ — k 2 = momentum sector under twisting in the g\ 
direction for parameters indicated by the black square in a). 

we confirm that the finite-size excitation gap of the KMS appears to remain finite in the 
thermodynamic limit. 

As we vary the magnitude of /5i and /5 2 , keeping all other parameters fixed, the spectral 
gap Ei — E changes dramatically (Fig. SI (a)), opening up as either amplitude approaches 
1 (corresponding to linearly polarized MW driving). At \/3\\ = \/3 2 \ = 1, diagonalization 
reveals a unique ground state, residing in the ki = k 2 = sector. While our system sizes are 
too small to clearly observe the Goldstone mode of a superfluid (SF), in contrast to the KMS, 
twisting the boundary conditions dramatically alters the ground state energy, as depicted 
in Fig. Sl(c); this is consistent with a SF which harbors long-range phase coherence. 



2 



